Draft version December 19, 2008 

Preprint typeset using WT^j^ style cmulatcapj v. 4/12/04 



RELATIVISTIC ACCRETION MEDIATED BY TURBULENT COMPTONIZATION 

Aristotle Socrates^'^ 

Draft version December 19, 2008 

ABSTRACT 

Black hole and neutron star accretion flows display unusually high levels of hard coronal emission 
in comparison to all other optically thick, gravitationally bound, turbulent astrophysical systems. 
Since these flows sit in deep relativistic gravitational potentials, their random bulk motions approach 
the speed of light, therefore allowing turbulent Comptonization to be an important effect. We show 
that the inevitable production of hard X-ray photons results from turbulent Comptonization in the 
limit where the turbulence is trans-sonic and the accretion power approaches the Eddington Limit. 
In this regime, the turbulent Compton ^/-parameter approaches unity and the turbulent Compton 
temperature is a significant fraction of the electron rest mass energy, in agreement with the observed 
phenomena. 

Subject headings: accretion - black holes, neutron stars 



1. GENERAL ISSUE AND MOTIVATION 

Like the envelopes of late- type dwarfs, accretion disks 
are thought to be turbulent and optically thick. Turbu- 
lence in stellar envelopes is driven by convection, while 
in accretion disks turbulence is commonly thought to be 
driven by the magnetorotational instability (MRI; Bal- 
bus & Hawley 1991). 

Accretion disks found in proto-stellar and proto- 
planetary systems, cataclysmic variables. X-ray binaries 
and quasars/ AGN all display coronal behavior as do late- 
type dwarfs. By "corona," we mean the appearance, in 
a spectral energy distribution, of a relatively energetic 
population of particles that persistently emits radiation 
at temperatures in excess of the given object's thermal 
photosphere. 

In the case of black hole and neutron star accretion 
flows, the ratio of coronal to bolometric power L^jL^^^ ~ 
0.1 — 1. In other systems that display radiatively effi- 
cient accretion, the value of L^/ L^^^ is smaller by orders 
of magnitude. For example, turbulent stellar envelopes 
acheive L^jL^^^ < 10^"^ in the most extreme (rapidly 
rotating) case. 

Here, we ask why black hole and neutron star accre- 
tion flows display such abnormally high levels of coronal 
power, observed in the form of hard power-law X-ray 
photons. 

In the directly observable case of the solar corona, me- 
chanical energy leaks out through the photosphere in 
some combination of waves and relatively steady mag- 
netic structures. The two reservoirs of stellar mechanical 
energy are convective turbulence and differential rota- 
tion. 

Since proto-planetary and white dwarf accretion disks 
as well as rapidly rotating late-type dwarfs display rel- 
atively modest levels of coronal power, differential ro- 
tation cannot be the sole key ingredient in explaining 
the extreme coronal phenomena seen in relativistic ac- 
cretion flows. Furthermore, the difference between MRI- 
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driven and convective turbulence must be ruled out as 
well because the non-relativistic accretion flows display 
only modest coronal phenomena. 

An obvious feature that is unique to black hole and 
neutron star accretion flows is that close to the central 
sources of gravity, where most of the accretion power 
is generated, the bulk motions approach relativistic val- 
ues. If the electrons possess bulk motions that are trans- 
relativistic, then they may readily amplify soft photons 
up to high energies, potentially forming a hard power- 
law X-ray spectra though the process of second order 
bulk Comptonization.'^ 

In the brief discussion that follows, we focus solely on 
the importance of turbulent Comptonization in relativis- 
tic accretion flows. Our approach differs from that of 
Socrates et al. (2004) and Thompson (2006) in that 
we completely ignore the relationship between the re- 
gion from which soft thermal seed photons originate - 
presumeably a cool thermal disk - with the region of 
Comptonization. We start with a description of an iso- 
lated eddy undergoing turbulent Comptonization. 

2. THE BASIC UNIT; AN EDDY IN A BOX 

First, we summarize the properties of an eddy in a box 
of dimension L. We assume that the eddy is embedded 
in an isotropic turbulent ffow. The scale of the eddy 
A < L is of order the box size and the characteristic 
eddy velocity v\<c\% taken to be trans-relativistic. We 
assume that the box is fully ionized with a constant mass 
density p. 

The eddy is stirred on the outer scale A ~ L and in 
the absence of any form of microphysical dissipation, the 
energy injected into the scale A flows down to smaller 
scales via a self-similar turbulent cascade. At the dissi- 
pation scale, the cascade cuts off and the turbulent power 
is converted into heat. 

We now discuss some of the relevant properties of our 
eddy. 

2.1. Some Important Properties of the Eddy 

^ As an aside, Cosmologists refer to this radiative mechanism as 
the "second-order kinetic S-Z effect." 
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The total kinetic energy 1C\ in the box at any given 
moment is given by 

Kx^ Mxvl^ pW'vl. (1) 

The total eddy power, or luminosity, L\ ^ IC/tx where 

tx - X/vx (2) 

is the eddy turnover time. We write the eddy luminosity 
as 

Lx^pvlX" (3) 

which is just the turbulent energy flux, multiplied by the 
eddy cross section. Also, we assume that in the absence 
of dissipation, the eddy velocity amplitude scales as a 
power law i.e., 

vx (A/Ao)" . (4) 

For isotropic incompressible Kolmogorov turbulence, the 
the index n = 1/3. We take vo < the speed of hght c and 

Since the velocity amplitude vx of our eddy is trans- 
relativistic, the inverse-Compton effect on electrons may 
become important on timescales of interest, particularly 
when the eddy is optically thin. Now, we assess the im- 
portance of the various microphysical processes arising 
from turbulent Comptonization. 

2.2. Eddy Comptonization: Cooling Rate and Compton 

Power 

The eddy cooling rate is the rate at which photons 
remove energy from the turbulent electrons, which are 
"stuck" to the ions as a result of - the assumed - tight 
Coulomb coupling. Therefore, inverse-Comptonization 
off of turbulent eddies may be viewed as a microscopic 
dissipation mechanism for the turbulence. Note however, 
that this dissipation does not necessarily occur on scales 
smaller than X ^ L. If the eddy optical depth ta < 1 
i.e.. 



spX<l 



(5) 



then the photons sample the turbulent velocities on the 
outer scale of the turbulent eddy, which is ^ vx ~ vq. 
Thus, the act of electrons, and thus the eddy, Compton 
up-scattering relatively cool photons necessarily implies 
that turbulent Comptonization off an eddy is a damping 
mechanism. 

Comptonization preserves the number of quanta, 
rather than the energy per quanta. If the radiation field 
is close to isotropic, the Kompaneets (1957) equation, 
which expresses conservation of photon number, reads 



3c2 



C6) 



where is the photon occupation number. The above 
expression is accurate as long as v\ 3> kBTe/nie and 
v\ ^ E^/rrif. = hv/rue- In other words, turbulent Comp- 
tonization is important as long as the contribution to the 
electron kinetic energy resulting from turbulence exceeds 
the contribution from random microscopic (e.g., thermal) 
motions as well as the average "seed" photon energy. 
Immediately, we can read off a time scale from eq. ([6]) 



where tc^x is the time it takes for the eddy to boost the 
energy density of the photon field by a factor of two. 
It follows that the turbulent Compton power of a given 
eddy is simply 



Lc,x - 



acsne c^UyX^. 



(8) 



If at any point along the cascade, the turbulent Comp- 
ton power Lc.x approaches the turbulent eddy power Lx, 
then turbulent Comptonization off of an eddy may be 
viewed as a viable microphysical dissipation mechanism 
for the turbulence. 

2.3. Electron Heating Rate 

There are really three particle distributions of concern 
in our box: the photons and the bulk turbulent and mi- 
croscopic thermal degrees of freedom for the electrons. 
So far, we have quantified the rate at which energy is 
transferred from the turbulence to the Comptonized ra- 
diation field. If the electrons are "cold" such that their 
thermal motions are small in comparison to their bulk 
motions, then the electrons stand a chance of heating 
up. 

In a frame moving with velocity vx , co-moving with a 
given portion of an eddy, the thermal energy density of 
the gas Ug evolves as'^ 



dUg 
dt 



aosUe C 



E-y 



rUeC^ 



(9) 



In the above expression, E^ is mean photon energy, av- 
eraged over the photon energy spectrum. The above ex- 
pression conveys the action of thermal Comptonization.^ 
We are interested in the regime where turbulent Comp- 
tonization boosts relatively few photons up to high en- 
ergies, where E^ < rUeC^- Though few in number, their 
large energies allow them to carry away a significant - if 
not dominant - fraction of the photon power. If the gas 
is cold such that the bulk - and therefore, turbulent - 
motions are large in comparison to the particle's thermal 
motions, then E^ ^ ksTe- In that case, we identify a 
Compton heating rate 



t 



E-y a 



c,hcat 
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(10) 



t: 



(7) 



which is the rate at which energetic Comptonized pho- 
tons heat up up the gas, increasing its thermal content 
per unit volume. 

3. AN EDDY AROUND A BLACK HOLE 

Now, we place our eddy near a black hole to get a feel 
as to how an ensemble of eddies - that collectively form 
an accretion flow - behaves. We place our eddy at a 
radius R = Rc/e from the hole, where Rg = GM,/c^ 
and e is the gravitational radius and radiative efficiency 

* Throughout, wc assume that the electron and proton temper- 
ature are equal to one another. 

^ One of the most famous examples of thermal Comptonization 
maintaining thermal equilibrium between radiation and gaseous 
matter is in the cosmic background radiation, before and during 
the moments of recombination. 
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of the hole, respectively. The radiative efficiency e is 
roughly determined by the location of inner-most stable 
circular orbit Risco i-e., e ~ Rg/Risco- So, if the e ~ 0.1 
then roughly ~ 100 MeV worth of gravitational binding 
energy per baryon is available for conversion into other 
forms of energy. The discussion that follows applies to 
neutron star accretion flows as well. Since black holes do 
not have a surface, the overall physical system is simpler 
and for this reason, we restrict our discussion to black 
hole accretion flows. 

3.1. Accretion power and accretion timescale 

We imagine that our eddy is surrounded by other ed- 
dies of comparable size and shape. The interaction be- 
tween eddies not only contributes to the turbulent cas- 
cade of energy to small scales, but also transfers angular 
momentum between radially adjacent eddies. 

The transfer of angular momentum, the viscous gener- 
ation of accretion power and the release of gravitational 
binding energy are connected by the expressions that 
govern angular momentum and energy balance. The vis- 
cous dissipation rate per unit volume resulting from 
an effective turbulent viscosity is given by 



£+ ^ PVI^^ pvln, (11) 



dlni? 



which is just the kinetic energy density, multiplied by the 
shear rate. Note that in the above expression, we assume 
that each respective eddy is rotating at the Keplerian 
rate, with an angular velocity fl ~ ^/GM^/R^ . The 
total accretion power, or accretion luminosity Lace, a for 
the eddy is given by 

Uoc,x^ P>?vln^JCxn. (12) 

Of course, each eddy is a transient structure embed- 
ded in a continuum of matter that slowly drifts towards 
the hole. This drift, which results from the transport of 
angular momentum is quantified by the specific angular 
momentum i and the mass accretion rate through the 
eddy, Mx, in the following way 

Mxic^X'pvl, (13) 
where we have ignored the contribution from the inner 
boundary. As an interesting aside, we see that M\i ~ 1C\ 
or, the rate of angular momentum loss within an eddy is 
equal to the kinetic energy, ICx, of the eddy itself. The 
combination of the expression for viscous dissipation, eq. 
(|12p . and angular momentum conservation, cq. p3|l . al- 
low us to connect the accretion power Lace. a to the ac- 
cretion rate M\ through the eddy i.e.. 



GMxM, 

R ' 



(14) 



analogous to the well known result. 

The characteristic time for an eddy to deplete its mass 
iaee,A IS obtained by re-writing conservation of angular 
momentum, eq. (|13|) . as 

Mx^^ -^R^n~Mxvl~1Cx 



acc.A 



(15) 



which is often written in terms of an anomalous viscosity 
and a Shakura & Sunyaev (1973) a-parameter. 



3.2. An Eddington- Limited Eddy 

Now, we determine under what conditions the vari- 
ous dynamical, thermal and radiative timcscalcs become 
competitive with one another. Throughout, we specify 
that the turbulent motions on the stirring scale A are 
trans-relativistic - a necessary condition for turbulent 
Comptonization to be a powerful effect. Another condi- 
tion, is that the optical depth must at least be of order 
unity, in which case of the rate of photon escape is equal 
to the rate of energy gain, yielding a relatively fiat (in 
V F,j) and powerful Comptonized spectrum. 

We take the limit where A ~ i? ~ ^g/^ ^-nd S '^l ^ 
ec^. That is, we approximate that the eddy is in dy- 
namic balance with the black hole's gravitational field at 
R/Rg ^ e^^ gravitational radii from the system's grav- 
itational center. In order to calculate the eddy optical 
depth, or alternatively the mean density, we further as- 
sume that the eddy is radiating at the Eddington limit 
for the eddy Lcdd.A i-c. 



GMxM, 

R 



e Ma 



GM,c 



^cdd,A 



.(16) 



Note that Ledd.A differs from the Eddington limit Ledd 
for the entire flow by a factor of 1/Att since the eddy 
on subtends ~ 1 steradian of the solid angle around the 
black hole. Putting together the expression above with 
eq. (jlSp we may write the eddy optical depth as 



TA 



-1/2 > I 



(17) 



At the Eddington limit, the Thomson optical depth 
across the eddy tx is of order unity and therefore, the 
photon field "samples" the energy bearing eddies on the 
outer scale A. Therefore, the rate for turbulent Conipton 
heating t~\ and the turbulent Compton luminosity Lc,x 
should be evaluated by using the outer scale length A and 
the outer scale velocity vx- Finally, for an Eddington- 
limited eddy, we find 

*A ^ - tc,\ - Oe,A - *dyn' (1^) 

where t^yn ^ Cg/L ~ 1/Sl. Furthermore, 

LcX — Laee.A (19) 

for an Eddington-limited eddy. 

The two expression, eqs. (fT8|) and (fT9|) above, are in- 
teresting results. Together, they indicate the release of 
gravitational energy resulting from turbulent viscous dis- 
sipation is equally matched by the Compton power that 
is generated by the turbulence as well. Therefore, it is 
possible that all of the fundamental physical mechanisms 
that mediate accretion i.e., the transport of angular mo- 
mentum, randomization of gravitational binding energy 
and the conversion of that energy into photon power is 
entirely determined by a single physical phenomena: tur- 
bulence. 

However, our picture is not complete. Bathing the 
eddy in such an intense (self-generated) hard radiation 
field has consequences for the thermal motions of the 
gas. Consider the form of the Compton heating rate, 
given by eq. (fTU]) . The act of turbulent Comptonization 



boosts the mean photon energy to a value E^ 
which implies that 



< 



c.hcat 



u. 

Ua 



7 +-1 



t7i > t: 



(20) 
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Fig. 1. — Cross section of critical accretion at the Eddington 
Limit L J, J J onto a black hole with radiative efficiency e = 1. The 
flow is divided into "Thomson units" in order to depict the fact 
that the number of protons A^p responsible for fueling an Edding- 
ton's worth of power in one light crossing time ~ i?G/c is given by 
the ratio of the area of the black hole to the Thomson cross section. 
That is, Np = A^e ~ Rq/t'^ where = e^/rrieC^ is the classical 
electron radius. In order to maintain a high radiative efficiency, 
e ~ 1, the orbit of each individual proton must undergo a deflec- 
tion of order ^ one Rq, which implies that Ari^/c^ ~ ^x/"^ ~ ^' 
Together, this implies that the turbulent Compton y-parameter 
yx ^ 1, leading to an observed flat spectral index F ~ 2. 



for an Eddington- limited eddy since U-y » Ug even if 
ksTg ^ Ey. This results from the fact that the radiation 
pressure is roughly a factor mpC^ /kBTp larger than the 
gas pressure, at the Eddington Hmit. 

The large value of the heating rate t^l^^^^^ leads to the 
conclusion that in equilibrium, the thermal energy per 
baryon equals the average energy of the hard Comp- 
tonizcd photons. In that case, the source term on the 
right hand side of eq. ^ approaches zero and thermal 
equilibrium for the gas is reached. 

Altogether, the overall picture for an Eddington- 
limited eddy is different from the viewpoint presented 
in Socrates et al. (2004). Apparently, it seems unlikely 
that gaseous species can persist at specific thermal en- 
ergies that are significantly smaller than w^. Therefore, 
turbulent Comptonization in black hole and neutron star 
accretion flows may also be viewed as an electron heating 
mechanism as well. 

4. THE TURBULENT COMPTON Y-PARAMETER AND 
THE MEANING OF ITS VALUE NEAR A BLACK HOLE 

We continue our assessment of the Eddington-limited 
eddy near a black hole introduced in the last section. 
The spectral features of saturated and mildly unsatu- 



rated Comptonization are well described by the slope of 
a self-similar Comptonized spectral feature and a high 
energy cutoff. The slope is often expressed in terms of 
the Compton y-parameter, given by the fractional en- 
ergy shift per scattering multiplied, by the average num- 
ber of scatterings. When y ~ 1, the "input" spectrum, 
responsible for providing the soft seeds, is deformed by 
the addition of a hard self-similar power-law component. 
In such an event, the ratio of hard Comptonized to soft 
seed power is also of order unity. When the eddies are 
optically thin such that < 1, we define the turbulent 
Compton y— parameter as 

The important result above, that the turbulent Compton 
y— parameter is close to unity, is not an accident. 

Consider the simple representation of extreme quasi- 
spherical black hole accretion given by Figure [TJ For 
illustrative purposes, the hypothetical case where accre- 
tion is 100% efficient i.e., e = 1, is depicted. Thus, dy- 
namical effects arising from general relativity, such as the 
presence of the last stable circular orbit, are completely 
ignored. Critical accretion, which implies t~^^ ~ ^dyn' 
while unit efficiency, e = 1, implies that X R ^ Rq. 
For this case, the number of protons Np responsible 
for fueling an Eddington's worth of accretion power is 
roughly given by Np ^ A,/aT ^ R^/r^, where A, and 
Te are the surface area of the black hole and the classical 
electron radius, respectively. 

It follows that Comptonized photons remove the accre- 
tion power of each individual proton on average scatter 
once with an electron - as long as the number of elec- 
trons equals the number of the protons. Or in other 
words, ta ~ 1 for highly efficient critical accretion. At 
the same time, the value of randomized proton velocity 
squared Ai;^ ~ ~ , since of order uip worth of 
binding energy is removed per proton. The combina- 
tion of these two attributes of critical efficient accretion 
onto a black hole yields the result that y^ ~ 1. Hence, 
Eddington-limited relativistic accretion is a natural tur- 
bulent Compton amplifier. 

5. SOME OPEN QUESTIONS AND CONCLUSIONS 

At radii immediately further out from the inner-most 
stable circular orbit, the ability for trans-sonic turbu- 
lence to Comptonize soft photons diminishes since 
decreases, while t\ increases. The latter implies that 
the input photons can only "sample" eddies of lower en- 
ergy, on smaller scales (Socrates et al. 2004). Also at 
larger radii, the ability for the flow to absorb a photon 
via free-free and bound-free processes increases (Socrates 
et al. 2006). We suspect that these adjacent outer re- 
gions serve as the source of - presumably thermal -soft 
input photons for the centrally located Comptonizing ed- 
dies. In the case accretion onto a neutron star, thermal 
emission from the stellar surface is additional source of 
soft seeds. 

We have yet to specify the actual source of trans-sonic 
turbulence itself. In the familiar language of the Shakura 
& Sunyaev (1973) thermal disk model, we have assumed 
throughout that en < 1 i.e., < c^. It is now widely 
accepted that the MRI (Balbus & Hawley 1991) gener- 
ically leads to a turbulent accretion stress, allowing for 
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the transfer of angular momentum to occur in highly 
conducting fluids. However, simulations of the MRI typ- 
ically find low values for the Shakura & Sunyaev (1973) 
a-parameter such that a ~ 0.1 — 0.01, though a clear un- 
derstanding for the value of the MRFs a is not currently 
in hand. 

Our main results for Eddington-limited trans-sonic 
turbulence near a relativistic source of gravity are as fol- 
lows: 

1) Turbulence is responsible for transporting angu- 
lar momentum while simultaneously converting gravita- 
tional power of protons into Comptonized photon power. 

2) The dynamical, accretion, Compton cooling and 
eddy overturn time are all comparable to one another. 

3) The turbulent Compton power of each individual 
eddy is comparable to the accretion power of that eddy. 

4) Electrons will inevitably heat up to the a signifi- 
cant fraction of the electron virial temperature. In other 
words, turbulent Comptonization is an important elec- 
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tron heating mechanism. 

5) The y— parameter due to turbulent Comptonization 
is close to unity. Therefore, the outgoing hard X-ray 
spectrum is a power law with flat spectral index F ~ 2, 
in agreement with observations. 

6) The results above are independent of black hole 
mass. 
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